home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 1998 November / IRIX 6.5.2 Base Documentation November 1998.img / usr / share / catman / p_man / cat3 / complib / CHIDI.z / CHIDI
Text File  |  1998-10-30  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. CCCCHHHHIIIIDDDDIIII((((3333FFFF))))                                                            CCCCHHHHIIIIDDDDIIII((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CHIDI   - CHIDI computes the determinant, inertia and inverse of a
  10.      complex Hermitian matrix using the factors from CHIFA.
  11.  
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  14.       SUBROUTINE CHIDI(A,LDA,N,KPVT,DET,INERT,WORK,JOB)
  15.  
  16. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  17.      On Entry
  18.  
  19.      AAAA COMPLEX(LDA,N)
  20.         the output from CHIFA.
  21.  
  22.      LLLLDDDDAAAA INTEGER
  23.         the leading dimension of the array A.
  24.  
  25.      NNNN INTEGER
  26.         the order of the matrix A.  KVPT    INTEGER(N)
  27.         the pivot vector from CHIFA.
  28.  
  29.      WWWWOOOORRRRKKKK COMPLEX(N)
  30.         work vector.  Contents destroyed.
  31.  
  32.      JJJJOOOOBBBB INTEGER
  33.         JOB has the decimal expansion  ABC  where
  34.         if  C .NE. 0, the inverse is computed,
  35.         if  B .NE. 0, the determinant is computed,
  36.         if  A .NE. 0, the inertia is computed.
  37.         For example, JOB = 111  gives all three.  On Return Variables not
  38.      requested by JOB are not used.
  39.  
  40.      AAAA contains the upper triangle of the inverse of
  41.         the original matrix.  The strict lower triangle
  42.         is never referenced.
  43.  
  44.      DDDDEEEETTTT REAL(2)
  45.         determinant of original matrix.
  46.         Determinant = DET(1) * 10.0**DET(2)
  47.         with 1.0 .LE. ABS(DET(1)) .LT. 10.0
  48.         or DET(1) = 0.0.
  49.  
  50.      IIIINNNNEEEERRRRTTTT INTEGER(3)
  51.         the inertia of the original matrix.
  52.         INERT(1)  =  number of positive eigenvalues.
  53.         INERT(2)  =  number of negative eigenvalues.
  54.         INERT(3)  =  number of zero eigenvalues.  Error Condition
  55.  
  56.      AAAA division by zero may occur if the inverse is requested and  CHICO  has
  57.      set RCOND .EQ. 0.0 or  CHIFA  has set  INFO .NE. 0 .  LINPACK.  This
  58.      version dated 08/14/78 .  James Bunch, Univ. Calif. San Diego, Argonne
  59.      Nat. Lab.  Subroutines and Functions BLAS CAXPY,CCOPY,CDOTC,CSWAP Fortran
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCHHHHIIIIDDDDIIII((((3333FFFF))))                                                            CCCCHHHHIIIIDDDDIIII((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      ABS,CABS,CMPLX,CONJG,IABS,MOD,REAL
  75.  
  76.  
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.